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FOREWORD 


Under Contract DAHC-04-69-C-0056 between IIT Research 
Institute (LITRI) and the U, S. Army Research Office-Durhan, 
LITRI is performing a study of fragmentation hazards to un- 
protected personnel. This is the final report of the study 
which was carried out during the period July 1971 to 
December 1971. 


The work was conducted for the Department of Defense 
Explosives Safety Board and supervised by Col. W. Cameron III, 
Chairman, Lt. Col. J, D. Coder, Project Manager, Dr. T. A. 
Zaker, Explosives Scientist and Mr. R. G. Perkins, Safety 
Engineer. Principal contributors to the presently reported 
work include L. A. C. Barbarek, D. I. Feinstein, H. H. Nagaoka, 
J. D. Rouse and J. R. Wingfield. Special acknowledgement is 
made to Dr. Zaker for providing trajectory algorithms and over- 
all technical direction to this study. 
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(P sGavnac® A computer program has been developed and documented tha 
generates the information necessary to establish minimum separation dis- 
tances between various munition types and personnel in order to mitigate 
fragment hazards. This information includes the fragment density at any 
point from an accidental detonation and the probability of injury to per; 
sonnel within the hazardous area, 


{' SUPPLEMENTARY NOTES 
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Published munition effectiveness data have been modified, utiliz{ 
ing statistical techniques, to emphasize a conservative approach, with 
respect to safety, by a better resolution of heavier fragments. The modi 
fied munition effectiveness data, which is input to the computer program 
have been included as a technical data appendix. 


' Computed results have been obtained for seven munitions and are 
presented in the form of contours of constant fragment density and 
damage probability in the horizontal plane surrounding the munition of 
interest. In addition, these contours have been simplified to yield an 
averse value of the hazard as a function of range from the accidental 
explosion. ~ 
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] ABSTRACT 
4 A computer program has been developed and documented 
ar that generates the information necessary to establish minimum 


separation distances between various munition types and per- 

a sonnel in order to mitigate fragment hazards. This information 
includes the fragment density at any point from an accidental 
detonation and the probability of injury to personnel within 
the hazardous area. 


Published munition effectiveness data have been modified, 
utilizing statistical techniques, to emphasize a conservative 
approach,with respect to safety, by a better resolution of 

= heavier fragments. The modified munition effectiveness data, 
which is input to the computer program, have been included as 
a technical data appendix. 


Computed results have been obtained for seven munitions 

we and are presented in the form of contours of constant fragment 
density and damage probability in the horizontal plane surround- 
ing the munition of interest. In addition, these contours have 
been simplified to yield an average value of the hazard as a 
function of range from the accidental explosion. 
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FRAGMENTATION HAZARDS TO UNPROTECTED PERSONNEL 
1. INTRODUCTION 


Existing quantity-distance standards for the manufacture, 
handling, and storage of munitions are based on the net weight 
of explosive filler contained in the devices in an unsubdi- 
vided magazine, or operating building unit. Safe distances 
are prescribed in tabular form essentially proportional to 
the cube root of the explosive weight. 


Because peak blast pressures from explosions of different 
yield are the same at distances scaled by the cube roct of 
the respective explosive weights, existing standards imply 
that the acceptable risk to a given target is based on a 
peak blast overpressure criterion alone. On the other hand, 
the field of fragments projected to the far field from accidental 
explosion of a munition store, consisting of inert munition 
component fragments and secondary fragments from any enclosure, 
does not satisfy the same similarity rules as does the airblast. 
Thus, defining an acceptable blast overpressure level at a 
target implies the acceptance of different levels of risk of 
damage by fragments, depending on the quantity and composition 
of the munition store and its enclosure. 


To develop quantity-distance standards based on con- 
sistent blast and fragment hazard levels requires determination 
of the damage risk due to fragments from accidental explosions 
as a function of the quantity and type of munitions, and of 
the characteristics of the source environment and the vulner- 
ability of the target. A previous report (Ref. rf was a 
first attempt aimed at applying engineering analysis, supple- 
mentary experimental efforts, and currently available data on 
fragmentation and damage criteria to the problem of estimating 
fragment hazards at explosives manufacturing and storage sites. 
While that report resulted in a fragment hazard model and a 
set of fragment density and damage probability maps for a 
variety of weapons and targets, this report is concerned with: 
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(1) Refining that model and making its corresponding 
computer algorithm usable by personnel moderately 
experienced at ballistic calculations, 


(2) Revising near-field fragment mass distributions 
to more accurately reflect their use as input in 
estimating fragment hazards, 


(3) Developing a new set of fragment density and 
damage probability maps; studying the effect 
of specific personnel protection criteria for 
unprotected personnel, and 


(4) Developing simplified two-dimensional relation- 
ships of fragment density and damage probability 
as a function of range. 


1.1 Problem Background 

Under Contract No. DAHC-04-69-C-0056 with the U.S. 
Army Research Office-Durham, LITRI has been conducting a 
series of investigations concerning fragment hazards associated 
with accidental detonation of munitions. This work has been 
performed under the direction of the Department of Defense 
Explosive Safety Board. 


Phase I of this study was concerned with establishing 
quantitative damage criteria in terms of fragment mass, 
velocity, and attack angle for various targets including 
standing personnel, vehicles, aircraft, buildings and open 
weapon stores. In Phase II an analytical model was developed 
to predict the density of fragments and the probability of 
damage to the targets considered in Phase I from explosion 
of individual munitions of various types. These included 
gun projectiles <... general-purpose bombs. Here damage 
probability contc.u.s were obtained in polar coordinates for 
a horizontal orientation of the munition axis in each case, 
Phase III attempted to extend the fragment hazard model for 
individual munitions to the case of multiple munitions in open 
stores (Ref, 2). The result was a limited demonstration that 
an analytic model could be developed to describe the initial 
fragment ficld of a stack of munitions. How?tver, it was 
also brought out that this initial fragment field was often 
related to munition case design, stack configuration and 
mode of initiation. 
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1,2 Program Objectives 


In the current research activity the intent has been 
to develop an analytic tool, usable by personnel moderately 
experienced at ballistic calculations, and capable of generating 
the information necessary to eStablish the minimum separation 
distance to personnel. The objectives of the study were: 


(1) To develop and document a computer algorithm 
which uses munition effectiveness tables as 
input and computes fragment densities and 
damage probabilities to unprotected personnel 
based upon the following criteria: 

e A hazardous fragment has a kinetic energy 
of 58 ft-lbs or greater, and 


e An acceptable density of hazardous fragments 
is not more than one per 600 sq ft. 


(2) To develop a rational scheme for revising published 
munition effectiveness data to more accurately 
reflect its use as input in estimating fragment 
hazards. 


(3) To develop a simplified means of relating 
fragment density and damage probability to a 
radial distance over a fixed sector of the 
ground plane. 


(4) To utilize the computer algorithm and the revised 
data in order to compute the fragment density 
and damage probability from the explosion of a 
single round of each of the following seven 
munitions: 
500 1b low drag bomb Mark 82 Mod 1 (H-6 load) 
750 1b Bomb M117A2 (Tritional Load) 
105mm Howitzer Shell Ml (Composition B Load) 
155mm Howitzer Shell M107 (Composition B Load) 
175mm Gun Shell M437A2 (Composition B Load) 
5"/38 Projectile Mark 49 (CCMP A-3 Load} 
8"/55 Projectile Mark 25 (Explosive D Load) 
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1.3 Program Accomplishments 


The major result of this study has been the de- 
velopment and documentation of a computer algorithm which 
will allow safety personnel, moderately experienced at 
ballistic calculations, to make the computations necessary 
to establish separation distances for personnel due to 
fragment hazards. A rational statistical scheme has also 
been developed and demonstrated which revises published 
munition effectiveness data to more conservatively reflect 
its use in estimating fragment hazards. A set of computa- 
tions were made for seven munitions and the results are 
presented as contour maps of total fragment number densities, 
damaging fragment number densities and injury probabilities. 
The contour maps have also. been simplified to yield curves 
of number density and injury probability as a function of 
radial distance from the munition source within a constant 
sector of the ground plane. (i.e., the nose, base and side- 
Spray sectors) These results serve both to demonstrate the 
capability of the computer algorithm and also as design aids 
for explosive safety personnel. 


1.4 Program Highlights 


The following sections of this report are organized 
in such a way as to first present the analysis on which 
the fragment hazard computer model is based, to next present 
an analysis of the munition effectiveness data which is input 
to that computer model and to finally summarize the study 
in the form of conclusions reached and recommendations as 
to model improvement. 


Appendices, covering the study's deliverable items, 
are included following the conclusion of the report. These 
appendices include a user manual for operation of the fragment 
hazard computer program, a listing of the computer program, 
the contours and simplified curves dascribing the fragment 
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seven munitions. The revised munition effectiveness data 
for the seven munitions were published as a separate classi- 
fied document. 


2 THE FRAGMENT HAZARD MODEL 


This section describes the analysis which forms the 
basis for the computational model for predicting fragment 
density and injury probability contours for various munitions 
The mathematical modei, illustrated in Fig. 1, is limited to 
the consideration of «he single munition without environmenta 
protection, The model has been programmed for digital computa 
tion, is modular and accepts as inputs: 


(1) The spatial distribution of fragment masses 
and velocities for singie munitions, which are 
defined for intervals of polar angle. (i.e., 
munition effectiveness tables) 


(2) The k-factors for the individual munitions, which 
express the relationship between fragment masses 
and projected areas for various munition types. 


(3) Vulnerability criteria for targets of interest. 
(e.g., kinetic energy and critical densities) 


Output of the model includes: 


(1) Fragment density contours showing distances 
to isodensity lines for all azimuths. Contours 
can be printed for all fragments or for various 
classes of fragments. 


(2) Injury/damage probability contours, showing 
ground distances to isoprobability curves at 
all azimuths for various munition/target 
combinations. 


Elements of this model include the following steps: 


(1) Munition performance data is converted to an 
interual form, All data is normalized into a 
a common internal form, wherein the average 
fragment velocity and the average mass and 
number density for each of several fragment 
mass categories are represented as tabular 
functions defined everywhere on the surface 
of s sphere of unit radius located at the 
origin, Gaps in the original date are filled 
in, and in addition the data may be smoothed 
if this seems desirable, 
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Fig. 1 FLOW CHART FRAGMENT HAZARD MODEL 


= (2) A terminal fragment field is computed on a 
a selected horizontal plane (i.e., the plane 
is above, below or passing through the. 
oe origin) utilizing routines for evaluating 
i fragment trajectories. The quantities 
«e computed for each mass category are impact 
velocity, impact angle and number of fragments 
per sq ft in a plane normal to the impact angle. 
These quantities are expressed as tabular 
functions of range and azimuth for output 
processing. 


ay 


me (3) Selected functions of the fragment field 

quantities computed in step 2 are computed and 

plotted in this step. A wide range of functions 
oe are available. The basic functions involve the 
fragment field alone, and do not consider charac- 
teristics of a target. Examples of these func- 
tions are number of fragments per sq ft, and 
total fragment kinetic energy per sq ft. 


The target functions use target characteristics 
to determine the number of damaging fragments 
and = e probability of damage at every point 

in the field. These functions use tables or 
formulas to determine whether or not a fragment 
is damaging, by finding the minimum velocity 
required by that fragment to damage the target. 
The fragment is considered to be damaging if its 
y velocity exceeds the threshold. 


° (4) Fragment field functions may be plotted either on 
the printer, or on an off-line plotter if the 
output tape is run through an appropriate post- 
processor. 


2.1 Wrajectory Analysis 


Large quantities of terminal ballistic property data 
are used in developing the outputs of the computational model 
described above. These data are generated from the equations 
i of motion for the fragments. Since these computations represent 
the bulk of the computational burden involved in exercising 
the model and their accurate evaluation is essential, it is a. 
desirable to utilize a highly efficient numerical procedure for 
calculating trajectories. 


Figure 2 illustrates the motion of a fragment moving 
‘ under the influence of aerodynamic drag and gravity forces in 


nonrotating local coordinates x, y tangent and normal to the 
trajectory. The corresponding equations cf motion are: 


Fig. 2 COORDINATE SYSTEMS AND TRAJECTORY 
GEOMETRY 


« 
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xX + pv x + g Sina = 0 (1) 


ytouv y +g cosa =0 (2) 


where dots denote differentiation with respect to tire t. In 
these equations, ¢ is the acceleration of gravity, v is the 
speed in the path, and a is the angle between the x - axis 
and the horizontal, Instantaneously we have x = v and y = 0. 


The aerodynamic coefficient 3 is given by 


1 = Cpw A/WN (3) | 
where Cy is the drag coefficient, w is the specific weight of | 
air, A is the cross-sectional area of the fragment normal to 
the rlight direction, and Wis the fragment weight. The fragment 
area and weight are related empirically (Ref. 3) through a ballis- 
tic density k as follows 

w=k abl? (4) 

t 
| 


In terms of k, the aerodynamic coefficient becomes 


b= Gp w/2 (20) 2/3 (5) 


An approximate local solution to the equations of mo- 
tion is obtained by separating the displacement into two parts, 
one a basic solution satisfying the local initial conditions 
and the equations of motion with gravity absent, and the other 
a pair of perturbations satisfying the linearized residual 


equations (Ref, 4). The results, applicable for small depar- 
tures of the trajectory from the local initial tangent, are 


equivalent to difference equations appropriate to an arbitrary 
time step in a numerical integration of the complete trajectory. 
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The displacement x is assumed to be of the form 


i 
q 
zd 
3 


x= XK, + a (6) 


where the basic solution X, satisfies 
x, + pxe = 0 (7) 


and the initial condition X =v, while the perturbation Xp 
satisfies the associated residual of Eq. (1). 


The drag coefficient C) in general cepends on the Mach 
number, and the atmospheric weight density w is a function 
of altitude, If both of these factors are assumed to be con- 
stant during the time interval of interest, however, the aero- 
dynamic coefficient ® is a constant and Eq. (7) is easily 
integrated, The results are 


x, = [log (1+u)] /P (8) 

x, = v,/(14) (9) 
where 

ust Pvgt (10) 


In these expressions, t is measured from the time at which the 
fragment is at the local coordinate origin in Fig. 1, and Vo 
is the value of v at that time. 


Substituting Eq. (6) and (7) into the equations of 
motion, expanding v in binomial series, and neglecting terms 
of second order and higher in i and y, we reach the following 
results: 


x. + 26K,x, +g sinu = 0 (11) 


y + x,y + g cosa = 0 (12) 
These equations are linear in the displacement perturbations 


x, and y, and can be integrated analytically by standard 
methods. The displacement and velocity perturbations are 


x, =7(e/2) t? sina (1+u/3)/(14u) (13) 
¥ =-(g/2) t? cosa[ u(i+u/2)-1log (14u) ] /u2 (14) 
%, = -gt sina [ 1t+u(14u/3) Y(14uy? (15) 
¥ = -gt cosa (1+u/2)/(1+u) (16) 
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The leading fac:r~s on the right in the foregoing equations 
express the position and velocity changes due to gravity in the 
elementary case of a drag-free trajectory. The multipliers con- 
taining u all approach unity as u vanishes, and can be viewed as 
corrections on the effect of gravity due to drag. 


] where u is defined as before by Ey. (10). 
? 


si The drag coefficient and atmospheric density are assumed 
to be constant during each time step at their values at the 
beginning of the step. The methol is self-starting in that the 
position and velocity changes are computed from initial values 
at the current step only. 


wer 


Initial values v = Vo and a= a, are assumed to be given 
at the fixed coordinate origin in Fig. 2. Equations (8), (9), 
and (13) through (16) give directly the displacement and 
velocity components after a typical time step t in the local co- 
ordinates, With respect to the fixed coordinates, the displace- 
ments during the time step are obtained from the relations 


( 
ax = x cosu - y sina (17) 
‘sy = x sina + ¥ cosu (18) 
while the rotation of the trajectory tangent is given by 
: -1 +/+ : 
fa = tan” (y/x) (19) | 
For low register trajectories, (i.e., those launched at 
angles less than that corresponding to maximum range) the above 
analytic perturbation equations with a = My furnish an approximate 
solution for the complete trajectory in one time step t, the 
total time of flight. At impact, the expressions for the position : 
; coordinates x and y are of the form: t 
x = (x, + xp) cosu,- ¥ sina, (20) 
.- y = (x, + xX) sina, + y cosa,= 0 (21) | 
= where x,, Xo and y are given by Eq. (8), (13), and (14), I 


2.2 Fragment Density Computation 


Munition effectiveness data are used to derive initial 
conditions for the ballistic trajectories of fragments, These 
data are in the form of initial velocity and number of fragments, 
in each of several mass categories, and are functions of polar 
angle measured from the nose of the munition. The resulting os 
fragment density at any point of interest can be computed 
deterministically from the known terminal points of fragments 


in all the mass intervals in each polar zone, 


An individual munition is regarded as a nonisotropic point 
source of fragments which is rotationally symmetric about its 
longitudinal (nose-to-base) axis. Thus, the properties of the 
fragments emitted by a single munition are functions of polar 
angle @ measured from the nose. The format of typical munition ef- 
fectiveness data (Ref. 5) is shown in Table 1. Fragments in all mass 
intervals are assumed to be emitted from a given polar zone at ; 
the same velocity. 


4 
{ 
i 
{ 
| 
{ 
{ 
j 


A single munition is assumed to be detonated at a level 
ground surface, with the munition axis horizontal. , In order to 
determine the probability of damage to targets at various dis- i 
tances and directions from the source, it is first necessary to 
calculate the number den. ties of fragments of different terminal 
ballistic properties in the field surrounding the source, 


The fragments from a munition, considered as a point 
source on the ground,may be regarded as originating from posi- 
tions on a hemispherical envelope enclosing the source. Each 
point of origin on the hemisphere defines an elevation angle | 
a, and an azimuth ¢ in the horizontal plane measured from the 
nose of the munition, 


In Fig. 3, let u, v, and w be rectangular coordinate axes 
centered at the source, with w the vertical direction and the 
nose of the munition in the positive v-direction. Let R, be the” 
(arbitrary) hemisphere radius, and @ be the pular angle measured 
from the nose of the munition. 
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Fig. 3 FRAGMENT SOURCE GEOMETRY 


Transformation formulas between the spherical coordinates 
(9,v) with the munition nose at the pole and the spherical co- 
ordinates (a,, é) with the pole vertical give the following 
relations: 


Oo 
} 


cos cos a, cos ¢ (22) 


cos ¥ sin a,/sin ) (23) 


Fragments on ballistic trajectories in the absence of wind 
travel in planes of constant azimuth ¢, 


For an individual munition, the initial fragment proper- 
ties are functions of polar angle @ only. Therefore, for given 
elevation angle a, and azimuth ¢, the associated polar angle 
can be computed from Eq. (22) and all initial fragment properties 
at that point on the hemisphere can be obtained from munition 
performance data, 


Now consider the fragments of a single mass interval, all 
having the same average ballistic properties, emitted from a 
particular munition at a point (a5) ¢) on the enclosing hemis- 
phere. For simplicity, consider target points in the ground 
olane of Fig. 3. The terminal point corresponding to this 
initial point on the hemisphere is uniquely determined for each 
mass category through the associated ballistic trajectory. 


The family of trajectories for fragments of a given average 
mass (i.e,, those in one mass interval) can be thought of as a 
mapping of points on the hemisphere into points on the ground 
plane. Fragments of one mass interval originating from an 
element of area Ro cosa Agia, on the hemisphere are projected 
into the element of area RA¢AR on the ground plane. Thus if 
ny is the number density of fragments in mass interval i at the 


Source hemisphere and n; is the number density at the terminal 


point R, we have 2 
R- cos a 
nN. = n° a OF 


i i 
R|dR/da,, | 
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This formula permits computation of the number density of 
fragments in each mass interval originating from a point on 
the source hemisphere at any azimuth ¢, since for given 
initial ballistic properties the terminal point R is known 
from trajectory calculations. To calculate the number density 
across a plane normul to the final trajectory, equation (24) is 
replaced by 7 
R? cos a, 
as oA (25) 
R|-sin Op da. 

where d¢ is the final trajectory elevation angle. 


For an individual munition the mapping derivative dR/da, 
is given by 


av 
dR_ _ _oR oR _g@ ) 
do. * Sa, * SV. Sa do (26) 
oO ) o o 
where Ve is the i.itial velocity. From Eq. (22) we have 
ag = sin a, cos ¢/sin 9 (27) 
Oo 


The partial derivatives 9R/ga, and dR/dV, can be obtained 
analytically for lower register fragments; however, for upper 
register fragments, they must be obtained by numerical differen- 
tiation of a precalculated ballistic file. The analytic ex- 
pressions for the lower register are 


aR L 2R - X cosa, A xy COS Og 2R-x,cosa,, es 
Dp ~~ tan Za, “sin 2a, tan a; 

a = - 2(R-x,, cosa, + Xo sina, /tana,)/V, (29) 
oO 


The derivative dv, /d@ is determined from munition effectiveness 
data, 


The computer model, moving along rays of azimuth in the 
horizontal plane, determines the terminal properties of lower 
register fragments at preselected radial increments of range 
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I out to maximum range. At maximum range the corresponding 
elevation angle is noted and the remaining elevation angle 

“ up to 90 degrees is divided into equal increments. The 

wh 


model, utilizing the multi-step trajectory routine with a 

ne fixed drag coefficient of 1.28, calculates the terminal prop- 
wr erties of the upper register fragments corresponding to each 
increment of elevation angle. The fragment number density can 
be determined exactly at prescribed terminal points for lower 
register fragments and, by interpolation, at these same points 
for the upper register fragments. These contributions are ac- 
cumulated to give the total fragment density at the preselec- 
gs ted radial distances. By symmetry it is only necessary to 

= perform the analysis in half the horizontal plane to one side 
- of the munition axis. 


oe 


aie For a particular target, nondamaging fragments based on 
Z the applicable target vulnerability criterion are excluded 
from the accumulated number density. The number densities 
so calculated represent, in a statistical sense, expected 
values of the numbers of impacts per unit area, That is, 
these results, calculated deterministically as a function of 
position in the horizontal plane, are the numbers of damaging 
ne impacts per unit area to be expected on the average. To 
calculate the associated probability of damage (i.e., impact 
by one or more damaging fragments) to a particular target 
requires a statistical representation of the process, incor- 


porating these expected values of fragment number density. 
In this study the target is a sphere of unit cross sectional | 
area and fragments are considered to intersect this target 


sa Peden 
ae 


across a plane normal to the final trajectory direction. > 4 


2.3 Target Damage Probability 


Calculation of the probability of damage to a particu- 
lar target requires a statistical representation of the process 
of impact by damaging fragments, involving the target area 
and the expected number of impacts per unit area. The number 
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density of fragments is a function of position of the 
spherical target in the horizontal plane of interest. 


f+ a 
e . L 


Consider a target of projected area A, normal to the 
trajectory of fragments of mass interval i at the associated 
terminal point. Let ny be the expected number of impacts per 
unit area of fragments satisfying the applicable target vulnera- 
bility criterion. The impact process is assumed to be uniformly 
random in the immediate neighborhood of the point of interest. 
That is, impact by a damaging fragment is equally likely on 
all equal elements of area in the vicinity of the point. 
Statistically, we then have what is termed a Poisson process. 
The probability p,(j) of exactly j impacts on the target by 


damaging fragments of mass interval i is given by 
= j = 

ice (n,A,) exp ( n,;A,;) 

PySJ 3! (30) 


The probability of impact by none of the damaging fragments of 
mass interval i is therefore 


py (0) = exp (-n,A,) (31) 


The target is assumed to be damaged if struck by one or 
more damaging fragments of any of the mass intervals. On this 
basis, the probability of damage to the target, q, is given by 


This last formula permits direct computation of the probability 
of damage to a given target at every point of the horizontal 
plane of interest from the accumulated number densities of 


damaging fragments and the associated projected areas. The target 


in this study was a standing man. The man was considered to 


present a projected area varying from 1,33 sq ft in plan to 9.0 sq 


ft frontal area, A, was then calculated as the sum of the pro- 
jections of both a vertical and a horizontal target on a plane 
normal to the final trajectory direction. ‘The vertical target 


had an area of 9.0 sq ft and the horizontal target had an area of 
1.33 sq ft. 
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2.4 Description of Model Output 


Contour maps, desc1.bing each of the seven wunitions of 
interest, listed in Section 1.2, were produced utilizing the 
computer model and are included as Appendix B, For each muni- 
tion there are three contour maps: one describing the fragment 
flux (i.e., target hits per sq ft) due to all fragments; another 
describing the fragment flux of all fragments with terminal 
energy exceeding 58 ft-lbs; and a third which presents the 
probability of serious injury to standing personnel computed as 
described in the previous secticn. 


Another set of figures is shown in Appendix C. Here 
again, there are three sets of figures per munition. These 
are simplified curves, derived from the corresponding contour 
maps, and represent fragment density and damage probability 
in certain sectors of the ground plane as functions of radial 
distance only. There are three of these fixed sectors: that 
is, 10-degree sectors in the base and nose of the munition 
and a third sector representing a peak side~spray. This third 
sector varied somewhat among all of the munitions, but was 
usually defined by a line of peak values appearing within 
+ 10 degrees of the 90-degree azimuth. This line of peak 
values did not usually include the outermost or lowest value 
contour, The average value of this contour in the side spray 
sector was determined by inspection from the maps. 


Botl. ©". cuntours and the corresponding curves represent 
the fragment fieid from a minimum range of 250 ft out to the 
prescribed limit of effect (i.e., one hit per 6000 sq ft for 
fragment densities and a probability value of 0.0001). The 
minimum range is a function of a prescribed input variable and 
was chosen as 250 ft in this study. Since the results represent 
a single unit of munition and it is anticipated that future use 


of these results may be directed at multipie unit stacks, it seems 


reasonable to assume that the far fields results will be of more 


interest and the minimum range of 250 ft will be quite sufficient, 
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3. MUNITION EFFECTIVENESS DATA he 


Previous studies (Ref. 2) have indicated published munition 
effectiveness data, in their present form, are not suitable i 
for far-field fragment hazard analysis. In general, these data 7 
place heavier fragments into one or two broad weight intervals i 
with a corresponding average weight. Since the heavier frag- yas 
ments travel greater distances and their resulting terminal bi 
kinetic energy effect is high, it is necessary that these frag- 
ments be resolved into an adequate number of weight intervals. 


r 
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This section of the report describes several tasks which 
were undertaken to provide a systematic approach to both illus- 
trating the need for and subsequently the methods utilized in 
revising munition effectiveness tables. 


Pan a OE” CONOEE GTO EN Te 


3.1 Documentation of Original Fragmentation Arena Data 


In collaboration with the Ballistics Research Laboratory : 
at Aberdeen Proving Grounds, the fragmentation arena test 
procedure and data analysis techniques used to obtain munition 
effectiveness data were reviewed. Arrangements were made to eee 
examine the original arena test firing records at the APG 
Technical Library for the following munitions: 


e 105 mm Howitzer Projectile M1 

o 155 mm Howitzer Projectile M107 

@ 175 mm Gun Projectile M437A2 

e 750 1b General Purpose Bomb M117A2 (APG Data) 


Subsequently, contacts were established at NWL to obtain similar 
arena test data for the following munitions. 


o 5-in/38 Projectile Mark 49 Mod O (VT) 

@ 8 in/55 Projectile Mark 25 Mod 1 (HC) 

e 500 1b Low Drag Bomb Mark 82 Mod 1 

e 750 1b General Purpose Bomb M117A2 (Eglin AFB Data) 
For each munition the following information was documented. 
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e Test munition physical measurements 

e Specifications for each arena test facility 

e Listing of mass groups and polar zones 

e Individual listing of fragment weights greater than 

150 grains. 

An understanding of the fragmentation arena test techniques 
is an important factor in the proper assessment of the larger and 
more hazardous fragments. The arena test is used to determine 
fragment mass, velocity and spatial distribution of high-explosive 
munitions, The munition axis of symmetry is located horizontally 
and is taken as the polar axis. Designating the polar angle by 
@ and the azimuth angle as ¢, the fragmentation characteristics 
are a function of 9, and gravity effects are assumed negligible. 
An arena is constructed of an appropriate size as illustrated 
in Fig. 4. The arena test is designed to sample fragmentation 
characteristics in various polar angle intervals ranging 
between @ = 0 degree at the nose, and @ = 180 degrees at the 
base of the munition. These data samples are used to predict 
the fragmentation characteristics for the entire munition. 


An important aspect of the arena test procedure is the 
rviationship between the sample size obtained from the incremental 
polar zone on the arena recovery panel and the corresponding 
munition polar zone. Calculations were made from actual arena 
setup data to estimate the magnification factor associated with 
the integrated munitions data, These data, summarized in Table 25 
indicate that only about 2 to 7 percent of the fragments are 
sampled in the 90 degree polar angle region (side spray) for each 
test. Thus, each fragment must be multiplied by a factor of 15 
to 60 to obtain integrated data. It should be nuted that the 
fragment sample size is about doubled at the 30 or 150 degree 
polar angle regions, and then increases rapidly to near 100 percent 
at 0 and 180 degrees (nose and base regions). 


Table 3 summarizes data associated with fragments in excess 
of 150 grains for the munitions of interest. The 150 grain 
fragment was assumed to be the minimum damaging fragment mass 
projected at low elevations. Table 3 indicates that fragments 
in excess of 150 grains represent only 5 to 18 percent of the 
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Fig. 4 FRAGMENTATION ARENA DIAGRAM 
(500 LB MK 82-1 BOMB) 
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groups beyond 150 grains are few in number and wide in weight 
interval, Since these heavier fragments are few in number and 
averaged over wide weight intervals, the data samples developed 
from arena tests are statistically inadequate in describing 

the characteristics for these potentially damaging fragments. 
As an example, the munition effectiveness data may list 60 
fragments in one polar zone, each weighing 1100 grains and with 
a common velocity, while all other zones are free from that 
category of fragments. As input to the hazard model, these data 
could yield a ring-shaped hazardous region in the far field. 

It may be more reasonable to assume these fragments to be dis- 
tributed statistically such as 30 unics at 900 grains, 20 

units at 1200 grains, and 10 units at 1500 grains, with the 
fragments being dispersed over three polar zones. In this 
regard, the distribution characteristics associated with frag- 
ments less than 150 grains could be used to predict the statis- 
tical distribution characteristics for the heavier fragments. 


It should be emphasized that the munitions effectiveness 
data are valid statistically for characterizing munition per- 
formance for fragments of interest to the user. In general, 
the heavier fragments are only considered in the arena test 
procedure to assure that a conservation of weight is maintained 
between the recovered fragments in the arena sample and the 
integrated munition effectiveness data. 


3.2 Revision of Munition Effectiveness Data 


Seven sets of data, corresponding to the seven munitions 
of interest, were considered and of these, six were altered 
Co overcome one of the following two fundamental deficiencies: 


1) Gaps - Toward the heavier mass categories, "data 
voids" appeared making it difficult to approximate 


the cumulative distribution of mass as a continuous 
function, 


2) Refinement - The cumulative distribution of mass 
loses definition where the final mass categories 
are an amalgam of all mass greater than the upper 
limit of the preceding category. 
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The first of these deficiencies were more predominate in 
the 175 mm and 155 mm munitions while the remaining sets of data, Pe 
with the exception of the 105mm, had poor refinement in the heaviest 
mass category. Each problem was handled differently, and in 
general, the problem of refinement was the most amenable since 
it was not necessary to add "fictitious" mass categories as it 
was in the case of data with void deficiencies. , 


It is prudent to comment that the best improvement which 
can be made on the experimental data is increasing the sample 
size. It is recognized then that with restricted data collec- 
tion techniques and Limited sample sizes, extrapolation and sta- 
tistical inference must be conservative and cautious. 


The basic indicators which were used to establish trends 


and to make decisions, were the average mass (m) per mass 
category, the associated standard deviation (0), the ratio 


-/m, and the average mass frequency (f) for each mass category. 


The primary data for a particular munition was processed 
tu print out a matrix of fragment weights and frequencies with 
mass categories ascolumns and polar zones as rows. Calculated 
inforraczion consisted of average mass m, over all polar zones j 
within each mass category i, the corresponding average frequency 
f,, the standard deviation ;, a ratio o,/m, for each mass 


category, the total weight W = =: Msp fige and the cumulative dis~ 
tribution of weight 
_k 2 m, sf, 
B= a jt Ke* 12-4540. N (33) 


for all N mass categories. The two fundamental deficiencies noted 
in the data were identified by inspecting the cumulative frequency 
curve for obvious gaps in the domain of the function and by 
observing anomalies in trend of standard deviation values for the 
mass categories. 


Two cases which will serve to illus rate the method, 


A. The 155 mm Munition 


The cumulative distribution curve (Fig. 5) was plotted and 
it was noted that although the average mass categories range up 
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to about 2800 grains, almost 100 percent of the weight has been 
accumulated at about 1200 grains. By inspection of the curve 
in Fig. 5 it was decided to start by adding a data point at 
about 2000 grains so that that portion of the curve would be 
less sparcely populated. The next decisions must phrase an . 
answer to the questions: how much, and in what manner? ao 


First and second moments were used as the indices of : 
central measure and dispersion of the distribution of mass ~ 
and frequency across the polar zones within a mass category. 

In general, the data available inthehigher mass categories 

is not sufficient to construct a reasonable frequency histo- 
gram. In the lower mass categories it becomes possible to do 
this and in order to obtain a graphical sense for the distribu- 
tion of this mass, histograms were constructed for several 
categories. Figure 6 depicts the weight distribution in the 
200-250 category for the 155 munition, There is another dis- 
tribution which will be used as a decision aid. This is the 
frequency of mass over all mass categories with polar zones 
constant. This curve is shown in Fig. 7. 


Essentially then, we have a two-dimensional polar zone - 
mass frequency matrix in which we would like to make additional 
entries, Toward this end trends are identified in two direc- 
tions: across columns (mass categories), and across rows 
(polar zones). The intersection of these trends at the point in 
question allows an estimate of the permissible. entries which 
can be made at this point. 


It was noted that the ratio 7, /a; tended to be fairly 
constant and, in case of the 155 mm munition, the average ratio 
was on the order of .061. Using this in conjunction with the 
intention of adding mass at ® = 2000, the dispersion measure is 
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estimated as .061 x 2000 = 122. The plot of polar zone fre- 
quency versus mass category (Fig. 7) was consulted, and a 
frequency of 30 was associated with the average mass of 2000 
grains, By observing the trend across the mass categories, 
attention focused on polar zone 7 (60-70 degrees) and polar zone 
11 (100-110 degrees). Further, by noting the distribution of 


frequencies across the polar zones in the previous categories, the 


same proportions are used to obtain: 


Mass Category m = 2000 grains 
Frequency f = 30 


Polar zones receiving entries = 7 and 11 


The frequency in previous \ 


it 
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mass category : 
m = 1230 | 

Distribution of mass in 

zones 7 and 11 for 

m = 1230 


ti 


79 


11/79 for zone 7 
68/79 for zone ll 


zones 7 and 11 for 


Distribution of mass in 
m ~ 2000 


= 11/79 x 30%4, zone 7 
68/79 x 30%26, zone 11 


— 


Now using o estimated at 122 for m = 2000, we obtain 


mto = 2122 
m- o = 1878 


It should be noted at this point that the selection of m, 7 

and conseciently m + * impose an implied condition of symmetry 
for the distribution of mass in that particular m category. 
This mass is further apportioned so as to agree with the trend 
generally defined by the preceding adjacent category thus 
bending the implied symmetry to conform with a distribution 
which is most Likely not symmetrical. When data is added in 
this way the recalculated average mass will not generaily azree 
with the original average mass m, which was a somewhat arbitrary 
selection to begin with. The new m will be fairly close and 

in case of the 155 munition the data point was added at 2090 
instead of m = 2000. The completed prefile for the new "fic- 
titious" mass category consists of the following: 


1. All polar zones empty except 7 and 11 
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2. A mass of 1878 grains with a frequency of 4 is 
entered in potar zone 7 


3. A mass of 2122 grains with a frequency of 26 is 
entered in polar zone 1l. 

To complete the addition of this new mass category we now 
consider the conservation of total mass effects. Obviously after 
the addition of a new mass the total mass T, has been increased 
by some increment Ty. Thus, we have from the original or 
primary data 


2 £555 7 Th (34) 


where Ms 4 is the fragment mass in the jth polar zone in the 


ith mass category, and fi. is the associated frequency. Mow 


with the added data points we have 
aatees aS * * S ws 


* * 
where ee and fh are the added data for the additional mass 


ates mule ene. 


categories m;. The objective is to scale dow all fie and 
i, proportionately such that for new faye 


Toward this end, let Pi; be that proportion of f.. which generate 


the new £54 such that Eq. (36) is satisfied for new fey 
Then for original fay: 
k N *k ot 


oa fs 4Ms 5S. p..f..m.. = : 
Gas ae ha rele Pig tiga; 1 (37) 


i=1 


Dividing both sides of (37) by Ty = 25 £i5Mi5> we have 


ae k..Mm.. =1 (38) 


recalling that x = 


the first term in (38) is p. Hence 
Ss P, fy am: 
= par ij ij ij 


aa oa =1 (39) 
Dek Leelee 
i={ “ij"ij 
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The result shows that we can readjust all the frequencies 
by one constant of proportionality. Hence, a complete set of new 
data is generated with an added data point on the cumulative curve 
with the conservation of mass. The new cumulative curve, plotted 
from the revised data set, is shown in Fig. 5. As in every other 
case, the curve is shifted downward, which introduces a conservative 
element into tne intci:pretation of the data, 


B. The 500 1b Bomb 


The cumulative distribution for this munition was plotted 
from the original data and is shown in Fig. 8. Without modification, 
the original data was divided into ten mass categories: 


0-10, 10-20, 20-40, 40-80, 80-120, 120-150, 
150-190, 190-230, 236-310, 310 + grains 


The problem in this case is poor resolution. The standard deviation 
for example is 3.26 for category number 4, 8.33 for category 

number 6, 7.93 for category number 8 and 698.5 for category number 10. 
Tr aidard deviation of 698.5 is not in itself unexpected, Because 
of «.5 magnitude, it is an anomaly since its ratio o/m = .754 

is tore than ten times the expected value of approximately .06. 
Moreover, it is a sharp break in the standard deviation trend, 

If we focus on mass category No. 10, then we find it labeled as 

310 +. This is misleading since, in fact, mass category number 10 
contains fragments which weigh as much as 5000 grains. To over- 
come this resolution deficiency, the 10th category was expanded into 
five additional categories such that the total number of mass 
categories became 15, 
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These mass categories were then added on to the data, replacing 
the former category number 10. The result was a better than 
50 percent improvement in the average o/m for categories 10 i 
thru 15 and a smoothing of the expected standard deviation trend, 
The altered cumulative distribution is shown in Fig. 8. 


The data received from the Navy was not directly usable with 
out extensive recoding and key punching. This was done for the 
750 1b bomb and the cumulative distribution appears in Fig. 9. 
This curve falls below the curves for the original data and its 
altered form, Since the lower curve implies more mass in the 
higher mass categories it also implies a greater damage influence 
range. In all cases the alteration of the data produced the same 
effect by lowering the cumulative distribution curve. 


The following is a summary of how the remaining data were altered. 


| No alterations were necessary for the 105 mm shell. 


750 1b bomb: originally had 12 categories. The label on 

category 12 was 450+. Six additional categories were added 

as follows: 

12 450-700 

13 700-900 

14 900-1200 

15 1200-1500 } 
16 1500-2000 


17 2000-3000 
18 3000 + 
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The cumulative distribution is shown in Fig. 9 
5 inch Munition: 11 mass 


categories were expanded to 14 as follows: 


11 250-450 
12 ° 450-700 
13 700-1500 


14 1500-2006 
The cumulative distribution is shown in Fig. 10. 
The 8 inch Munition: 10 mass 


categories were expanded to 16 as follows: 
10 160-300 

11 300-500 

LZ 500-700 

13 700-1000 

14 1000-2000 

15 2000-3000 

16 =©3000+ 


The 175 mm munition: two data points were added 

to the cumulative distribution as shown in Fig. 12. 
The technique was the same as described for the 155 mm 
munition. 


A complete set of the revised munition effectiveness tables 
are presented in Appendix E to this report. Their format is 
similar to that shown in Table l. 


3.3 Model Sensitivity to Input Data 


In order to gain an appreciation for the sensitivity of 
the fragment hazard model output to alteration of munition effective- 
ness data, which is input to the model, a computational experiment 
was conducted. The unaltered and altered data sets for the 175 mm 
gun shell were used as input to the model. The results of these 
two cases are shown respectively in Figs. 13 and 14, The primary } 
difference between the two sets of contours is the location of 


The cumulative distribution is shown in Fig. ll. 
| 
i 
i 
| 
| 
| 


the lowest value contour Line. 
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Fig. 13 175 MM GUN SHELL M437A2 (COMP. B LOAD) 
(UNALTERED DATA) 
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CONTOURS OF TOTAL FRAGMENT 
R DENSITY 


300 
Range, (ft) 


CONTOURS OF FRAGMENT NUMBER DENSITY 
FOR FRAGMENTS WITH FNERGY 
GREATER THAN 58 FT LAS 


1500 1000 500 U S00 1000 1500 


Rauve, (ft) 
lesa Tail 


CONTCAUKS OF PROBABILITY OF INJURY 
TO PERSONNES. 


1500 


Range, (ft) 
Tesh 


Fig, 14 175 MM GUN SHELL M437A2 (COMP, B LOAD) 
(ALTERED DATA) 
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The altered data set represents an emphasis on the resolution 
of the heavier fra- mts, (e.g., note Fig. 12) at the expense 
of lower weight fragments, This has the pronounced effect of 
shifting the location of the lowest value contour line outward 
for the altered data set. In terms of safety, then, the altered 
data set results in a more conservative set of contours, 


4. RESULTS, CONCLUSIONS, AND RECOMMENDATIONS. 
This study has resulted in the development of a computer 
model which generates the information necessary in establishing 
minimum separation distances between various munition types 
and personnel in order to mitigate fragment hazards. The 
model specifically treats the fragment hazard associated with 
a single munition and has been utilized to generate single unit 
fragment hazard data for seven common military munitions. While 
single unit detonation does not represent a realistically severe 
accident situation, previous work indicates that multiple unit 
(i.e., stacks) fragment hazards may be proportional to single 
unit results, If this is true, and there is evidence to support 
such a hypothesis in the case of thin wall munitions, then the 
results obtained herein are directly applicable to estimating 
the fragment hazards associated with openly stored munitions. 


4.1 Results 


e Contour maps, expressing number density and injury 
probabilities (i.e., Appendix B), are seen to be 
particularly pronounced in three fixed sectors of 
the ground plane; that is, the nose, tail and side- 
Spray sectors. Average values within each of these 
sectors have been expressed as functions of the radial 
distance from the explosive source, (i.e., Appendix C) 


e Examination of the relationships in Appendix C indicate 
that the side spray sector is the critical area in the 
case of thick wall "shell type" munitions and to a 
limited degree in the thin wall "bomb" munitions. 


e The curves in Appendix C also serve to illustrate how 
the kinetic energy criterion for personnel injury 
(i.e., 58 ft-1bs) reduces the applicable number 
density, At a given range the number density is re- 
duced as much 4s an order of magnitude. It shuuld be 
recognized that consideration of a less severe kinetic 
energy criterion will lead to lesser reduction in num- 
ber density and this a more conservative injury cri- 
terion in terms of safety. 


e Table 4 summarizes the results available in Appendix C 
for each of the seven munitions considered. The table 
gives the computed distances corresponding to a fragment 
density of one hit in 600 ft? for ail fragments and frag- 
ments with a terminal energy in excess of 58 ft lbs. 
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TABLE 4 


SUMMARY TABLE FOR THREE PRINCIPAL DIRECTIONS AT 1 HIT 
PER 600 SQ FT (ENTRIES ARE IN RADIAL FEET) 


er RS ye RRS TS 


All Fragments Hazardous Fragments sa 


Munition 

Nose Side Tail Nose Side Tail 

750 (1)* 440 1060 740 220 690 500 
500 (1) 220 825 595 210 670 450 
175 (1) 250 840 575 250 450 200 
155 (1) 290 810 510 120 400 230 
105 (0) 240 650 360 100 270 150 
8 (1) 325 660 240 140 520 120 

5 (1) 310 720 340 140 275 150 


*(0) = Unaltered Data 
(1) = Altered Data 


** 
Hazardous fragments are defined as having in excess of 
58 ft lbs kinetic energy 
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4.2 Conclusions 


* 


e By utilizing trajectory analysis in conjunction with 
stochastic treatment of experimental input data and 
damage functions, a mathematical model has been de- 
veloped for estimating injury/damage contours for 
various single unit munitions. 


e@ The mathematical model currently utilizes published 
munition effectiveness data which has been altered 
to give more emphasis to heavier fragments, It has 
been observed that this results in extending contour 
levels outward giving a more conservative result 
insofar as fragment safety is concerned. 


4.3 Recommendations 


Before making specific recommendations concerning the 
results generated in this study and future research efforts it 
is prudent to discuss the adequacy of munition effectiveness data 
as input to the model, 


The results generated by the computer model are dependent 
upor) ind quite sensitive to the munition effectiveness data, 
which is input. This data was originally generated to support 
munition effectivcuess studies and is the result of explosive 
tests of single unit munitions. It is the only known source of 
information concerning near-field estimates of munition fragment 
size, number and initial velocity. However, since it has been 
collected to be utilized in weapon effectiveness studies, it is 
primarily concerned with the fragments which are effective within 
the applicable range of the munitim. This has ‘normally led to 
a set of data which has a high degree of resolution, in terms of 
weight intervals, where the greatest number of fragments are con- 
centrated. This unfortunately is at a rather low fragment weight 
(e.g., below 300 grains). The remaining fragment weight of the 
munition is quite substantial, but because it does not break 
down into very many fragments and is not always projected into 
the designed zones of munition effectiveness, its recorded 
resolution is usually quite poor. 


Another inadequacy of recorded munition effectiveness data 
is concerned with its use in representing the basis for multiple 
unit munitions fragment hazard analysis. Here, the primary 
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concern is whether the munition fragment size, number and initial 
velocities will be similar for munitions in single and multiple 
units. 


In the present study an attempt has been made to overcome 
the first of these deficiencies in a conservative manner. That 
ig, existing munition effectiveness data have been altered to 
adequately resolve higher weight fragments in a statistical 
sense. Recently reduced data, developed at the Naval Weapons 
uaboratory, tend to support these alterations; however, the 
best way to resolve this problem would be the design and use 


of an experimental procedure aimed at obtaining arena data for 
hazard analysis. 


em, hc oN eT Ba 


act 


the problem in assuming similar fragment characteristics 
for multiple and single unit stacks is a much more difficult 
problem to resolve. A number of detonation tests of stacked 
munitions have been conducted in the past where resulting 
fragments have been collected. In some cases the fragments have 
also been sized and number and weight distributions computed, ac 
Results of these studies indicate that thinwall "bomb type" ii 
munitions tend to fragment into similar size fragments for both 
multiple and single units. H-wever, it is becoming quite apparent 
that this is not the case for thickwall "shell type" munitions. 
Here, fragment weight and number distributions from multiple 
units are quite different from single unit munitions, 


coer 


In light of the input deficiencies enumerated above, the 
following recommendations are made: oe 


@ To take care in utilizing results obtained in this 
study for single unit munitions, in projecting 
fragment hazards associated with multiple units of 
similar munitions, This is especially true of 
"shell type" munitions, 


@ To compare the analytic results for the single unit 
750 1b bomb obtained in the present study with experi- 
mental multiple unit results obtained in the NWc- 
China Lake Tests of March, 1970. Such a comparison 
will serve to validate our analytic procedures and 
the extension of our results for estimating multiple 
unit "bomb" fragment hazards, 
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e To characterize the resulting fragment number and 
weight distributions obtained from the June, 1970 
155 mm shell tests at Yuma and the corresponding a 
155 mm shell tests at China Lake in December of 
1971. These tests are expected to be dissimilar 
amongst both themselves and as compared to single 
unit results. However, such comparison will x 
document these differences and the resulting dis- 
tributions can be utilized by the enalytic model 
to estimate the fragment hazard associated with 
both these cases. 
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APPENDIX A 


PROGRAM USER'S MANUAL 


APPENDIX A 
PROGRAM USER'S MANUAL 


The computer program for the fragment hazard model has 
been written in such a manner as to accept a problem-oriented 
input language. The following discussion describes, in detail, 
the various control and data statements required to execute 
the program. An example of the use of these statements in a 
sample problem is included, 


Control Cards 


All control cards start with column 1 = '$', The remainder of 
the card is free format. Only columns 1-72 may be used for con- 
trol text; columns 73-80 are not examined and may be used for 
anything. 


Control Statement Form 
Control statements are of the form: 
name fieldl field2 ... fieldn 


where 'name' identifies the control statement, and the 'field's 
are parameters whose form depend on the particular control state- 


ment. name' and the 'field's are separated from one another by 


1) one or more blanks, or 
2) a comma, padded on either side by one or more blanks. 
Examples: SPRINT CON,DATA 
STAPE 4,SRCH=01100,BKSP 
SSTOP 


If desired, comments preceded by a 'S' may appear after the last 
field. 
Example: $STOP $ END OF RUN 


If a control statement will not all fit on one card, it may be 
continued to additional cards in the following way: 


1) place a comma after the last field on the first card. 
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2) 
3) 


4) 
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start the second card with a '$' anda _',' (padding 
with blanks if desired). 


place remaining fields after the '$' and ', 
second card. 


as many continuation cards as necessary may be used. 


' on the 


Example: SORDNANCE FORMAT=2 ,REDUCE,REZONE ,NORM,SMOOTH , 
$ ,READ=3 , LIST $ READ FROM TAPE 3 


Parameter Form 


Control statement parameters may be of three forms: 


1) 


2) 


3) 


Positional: A value which must appear in its proper 
order on a control statement, as the first parameter, 
third parameter, etc. 


Examples: 3, 14.197, T, ‘TITLE A' 


Flag: A name, or the two characters 'NO' prefixed to 

a name. A flag may appear anywhere in a control state- 
ment after the positional parameters (if any). A flag 
is associated with a logical value: this value is TRUE 
if the flag's name is specified, or FALSE if 'NO' plus 

the flag's name is specified. 


Examples: LIST, NOLIST, BCD, BIN 


Keyword: A name, followed by either a BCD (3-8 punch) 
or an EBCDIC (6-8 punch) equal sign (=), which may be 
padded with blanks, followed by a value. 


A keyword parameter may appear anywhere in a control 
Statement after the positional parameters (if any). 


The variable represented by the name is set to the pro- 
vided value. : 


Examples; ORDER=2, Z=5.5, FLAGS=XY, TITLE='TITLE A' 


Value Form 


There are four types of values: 


1) 


Logical: a string of one or more characters, not con- 
taining any of the following: dollar sign ($), com- 

ma (,), equal sign (=), blank, or apostrophe ('), If 

a T appears in the string before an F appears, the 
value is TRUE. If an F appears before a T appears, or 
if neither an F nor a T appear, then the value is false, 


Examples: T, .TRUE., TRUFFLE and TRUE are all TRUE. 
F, .FALSE., AFTER, FALSE and NO ere all FALSE. 
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2) Integer: one or more digits, optionally preceded by 
a + or - sign. 


Examples: 0, -1234, +77111 


3) Real: a real numbex. Decimal point and exponent may 
be used. The exponent, if present, must be of one of 
the following forms: 

En, Etn, E-n, +n, -n 
where n is one or two digits. 


Examples: 0, ©1234, +77111, 3.1416, 6.0238E+23, 1.-5, 
-.1+10, 51E6 


4) Alphabetic: two forms are allowed. 


@ a string of one or more characters, not containing 
any of the following: dollar sign "($) 3 comma (,), 
equal sign (=), blank or apostrophe ('). 


Examples; ABC, 123, MACH-1 


ea string of zero or more arbitrary characters, en- 
closed in apostrophes ('). Within the string, one 
apostrophe is represented by a pair of epost opucks 
mituer tue sCD (4-3 punch) or EBCDIC (5-8 punch) 
apostrophe may be used. 


Examples: 'ARTHUR', ‘WHICH WAY?', 'DON''T GO", 
"TITLE A’, ‘' ( a null string) 


Notation 


When describing control cards, some possibly unfamiliar notation 
is used. This can be explained with an example: 


S$TAPE “y aah ,REW] 


The use of { }means that only one of the arguments listed may be 
used. If one of the arguments is underlined, that one will be 
assumed if none of the group are specified. (In the example, BIN 
is assumed unless BCD is used. BIN and BCD may not both be used.) 
The use of [{ ] means that the enclosed argument is not required; 
it may be used or not used, 
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SPRIN1 _ ,CON][ ,DATA] 
Print controls 

CON Cause control cards to be printed. 

DATA Cause data cards to be printed 
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SSTOP 
End of run; execution is terminated. 


P-int elavsed time since beginning of execution, 
and since last STIME statement. 


BCD). a SRCH=a 
STAPE n (BIN ,REW]{ ,SKIP=k] {° xSRatnat ,BKSP}| ,WEOF] 
Tape file manipulations. If more than one of 

the operations (REW through WEOF) are specified, 
the operations take place in the order REW,SKIP, 
SRCH or SXRCH, BKSP,WEOF, regardless of the order 
in which the operations are specified. 


n Tape unit number 
BCD Tape is formatted (BCD) 
BIN Tape is unformatted (binary). 
If neither BCD nor BIN is specified, BIN is assumed. 
REW Rewind unit n. ; 


i $TIME 
| SKIP=k Skip k records on unit n. 
SRCH=@ Search unit n, starting at its present position, 
for a sentinel record whose key is a (1 to 8 
characters). Terminate if no such sentinel is 
found before the end-of-file. 


XSRCH=a Proceed as with SRCH, but if an end-of-file is 


ioe os 


| encountered, rewind unit n and search the entire 

file. Terminate if end-of-file is again reached 
without finding the proper sentinel. 
| BKSP Backspace unit n one record. 
WEOF Write a terminal sentinal record and an end-of- 
| file mark on unit n. 
SORDNANCE FORMAT=n{ ,NSETS=n][ ,REDUCE][ ,REZONE][{ ,NORM][ ,SMOOTH] 

| [ ,ORDER=n]{ ,READ=n]{ ,LIST][ ,KWIKPLOT]({ ,SUMMARY] , 
{ 


Cause a title card and a set of ordnance data 
to be read. 


FORMAT=n Ordnance format code. =1 or 2. 


NSETS=n Number of ordnance components (fuz_, casing, etc.) 
If<0 or omitted, NSETS=1 is assumed, 
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REDUCE Halve the number of mass categories by combining 
each pair of categories. 


REZONE Transform data from N+l sets of values at the 
boundaries of N polar zones to N sets of values 
at the N polar zone midpoints. 


NORM Normalize data, given as total number of frag- 
ments per mass category per polar zone, to 
number per unit area on a unit sphere surrounding 
the ordnance. 


SMOOTH Fill in gaps in data, and (if ORDER>O) apply a 
smoothing function to the ordnance data. 

ORDER=n Use an nth order Fourier series to smooth data 
(if SMOOTH is specified). 


READ=n Read ordnance data from unit n, a BCD card image 
unit. If n=5 or if READ=n is omitted, ordnance 
data is ussumed to follow the SORDNANCE statement. 


LIST List ordnance data in tabular form. s 


SUMMARY List data summary information: total mass per 
mass category, per polar zone, etc. 


KWIKPLOT Produce printer plots of velocity curve, and of 
mass and number curves for each mass category, 
all as a function of polar angle. 

An ordnance title card is read immediately after a SORDNANCE state- 
ment is encountered. Columns 1-72 contain a description of the 
ordnance. 


Two basic formats are defined for ordnance data (FORMAT=1 and 2). 
Within the context of each format, several data conditioning op- 
tions are available (REDUCE, REZONE, NORM, SMOOTH, ORDER=n). 


Ordnance data is read from cards if the READ field is omitted, or 
if READ=5. If a READ field specifies some other I/O unit, input 
is read from that unit in card image form, This unit is assumed 
to be properly positioned. 


The first card (or card image) of a set of ordnance data is as fol- 
lows, regardless of the type of ordnance data. 
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COLS _FORMAT NAME DESCRIPTION aes ae 
1-5 I5 NMORD Number of mass categories. 


T£ NMORD> 30 (maximum number of mass 
categories), REDUCE must be specified 
to cut the number of mass categories in 
half. 


6-10 I5 NVORD Number of measurement positions along a 
meridian from nose to tail. Measurements 
may be made either at the center or at 
the edges of equi-spaced polar zones. In 
the latter case, REZONE must be specified 


NVORD¢ 37. 
ll-lz 45 IDORD Ordnance {£.D. uumber O<iDORE<99. 
16 (ignored) 
17 Al Ivo Ordnance modification code. If left 
blank, IVO='0' is assumed. 
18-20 (ignored) 
21-30 £10.0 XKORD Average fragment drag factor (erains/in.> 


Format 1 Ordnance Data 


Data of format 1 consists of NVORD groups of cards, each containing: 


@ [NMORD/4]" cards, in (4(2F5.1,5X))format, containing NMORD 
pairs of values; the average Exgement mass (grains) and 
the average number of fragments’ Tor each mass category, 
for one polar zone. 


e one card, in (50X,E10.1) format, containing the average 
fragment initial velocity (fps) for one polar zone. 
Format 2 Ordnance Data 


Data of format 2 is split into two parts. 


Part 1 contains NSETS of cards, each consisting of [{NMORD/4] groups 
of NVORD cards. Each card is in (8E10.0) format and contains 
four pairs of values: the average fragment mass (grains) and 

the average number of fragments for four mass categories, for one 
polar zone. 


aie Neaiis che smallest integer not less than X. 
Examples: [3.14]=4, [14]=14, [7.00001]=/[7.9999]=8. 


Yragment number information in the form of a count of the tctal 
number of fragments in a polar zone may be used, provided that 
the NORM parameter is specified on the SORDNANCE statement. 
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Zone 1 A Zone 4 
Zone 2 Zone 3 
Hypothetical Sphere of Unit 


Radius Surrounding Ordnance. 


Sphere divided into four polar 
zones, each of equal latitude 


range. 
a yeaa 
; x 
: 
4 
i 
/ 
Y a“ 
A «eho a 
Measurements Taken at Center 
of Zones i 


>, : a 


Measurements Taken at Zone 
Boundaries (REZONE 


Parameter Required on 
SORDNANCE Statement) 
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Part 2 contains NSETS sets of cards, each consisting of NVORD cards 
J in (50X,E10.0) format, each of which contains the average fragment 
initial velocity (tps) for one polar zone. 


$BOOM | ,FLAGS=a]| ,NBAR=n]( ,ORDER=n}| ,Z=h][ ,DRNG=r][ ,NRNG=n] 
' NEL=n}[ ,NAX=n] 


Generate a fragment distribution using ordnance ig 
data read according to the prior SORDNANCE state- ; 
ment. Write the fragment distribution field, : 
with a sentinel record on binary unit 4. i 


AeA AR AR OS aR 
ms 


ae bh 


7 FLAGS=a « consists of one or two characters. If a is one 

J character, the second flag character is assumed : 

: to be zero. If the FLAGS field is omitted, both ‘ 
flag characters are assumed to be zero. The flag : 

if characters are used in the sentinel I.D. for the 

ey fragment field output. 


NBAR=n Number of barriers placed in the vicinity of the 
: ordnance. If omitted, NBAR=0 is assumed. Barrier 
description cards, if any, are read after the frag- 
ment field title card, NBAR<1. 


it “z=0 height of target plane (in feet) above ordnance. 
If Z<€0, the ordnance is above the target plane. F 
i DRNG=r spacing of range points, in feet. 
- NRNG=n number of range points on an azimuth ray, NRNG <20. 9 
’ NEL=n. Number of elevation angles to be considered (high § 
j register fragments), NEL<18. : 
NAZ=n number of azimuth rays in fragment field (a semi- 


circle), NAZ< 36. 


The following cards are read immediately after a $ BOOM statement 
is encountered: 


@ A title card. Columns 1-72 contain a description 
of the fragment field. 


e If NBARD)O, NBAR barrier description cards are read, 
The format of each card is: 
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COLS FORMAT NAME 
1-5 T5 I3TY?P 
5-10 

11-20 E10.0 BCPH 

21-30 E10.0 BWPH 

31-40 E10.0 BHAL 
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Circular Barrier 


gS a ee 


{i 
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DES CRIPTION 


Barrier type code. 

l=barrier is a circular arc, with 
ordnance at center. 

2=barrier is straight and is closest 
to the ordnance at its midpoint. 


(ignored) 

Azimuth angle (deg) of barrier 
midpoint 

Extent of barrier in azimuth (deg) 
on either side of midpoint. 


Angular height of barrier (deg), as 
seen from the ordnance. If IBTYP=2, 
this angular height is measured at 
the midpoint of the barrier. 


vv” 
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Straight Barrier 
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, TABLES 
»PARAMS 
;LIMITS 


»SOJAC 


;EXEC[ ,PLOT][ ,KWIKPLOT] { [ ,LIST] 
RANGE 


TABLES 


PARAMS 


LIMITS 


so 
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Read function definitions or compute 

and plot functions of the fragment ; 
field variables (velocity, striking i 
angle, and number per square foot, as : 
functions of range and azimuth). i 


The fragment field is read from binary 
unit 4, which is presumed to be prop- 
erly positioned. 


Read Damage Threshold Velocity Tables 
or Formula Coefficients. 


Between two (2) and five (5) cards 
are read for each Table or Formula. 
The formats of these cards are given 
below, The end of input is marked by 
a card containing blanks or zeroes in 
columns 1 and 2. 


Read a titie card, and fragment thresh- 
old and target description cards for 
each fragment field function, 


Between two (2) and four(4) cards are 
read for each function. The formats 

of these cards are given below. The 

end of input is marked by a card con- 
agri blanks or zeroes in columns 1 
and 2. 


Read an output parameter card for each 
fragment field function, 


The format of this card is given below. 

The end of input is marked by a card 

oe blanks or zerces in columns 
and 2. 


~~ 
The formula used is the empirical relationship derived under 
project THOR. 


EXEC Compute and plot specified functions 
of a previously computed fragment 
field. Functions are specified on the 
data card following this control card. 
The format of this card is given below. 


The fragment field is read from binary 
unit 4, which is presumed to be properly 
positioned, 


The functions are specified using two 
digit codes, whose meanings are: 


01 total number of fragments (per sq ft) 
02 total fragment mass (lb) 
03 total fragment momentum (lb-sec/ft 


04 total fragment kinetic energy 
(f£t-1b/f£t-) 


05-07 (unused) 


08" average fragment impact velocity 
(ft/sec) 


o9* average fragment impact angle 
(radians) 


lj number of damaging fragments, 
using table j. 


2j probability of damage, using table j. 


3j compute and plot 1j and 2j together 
(this is cheaper than specifying 1j 
and 2j separately) 


4j number of damaging fragments, using 
Thor coefficient set j. 


5j probability of damage, using Thor 
coefficient set j. 


6j compute plot 4} and 5j together 
(this is cheaper than specifying 
4j and 5j separately) 


7j-9j (unused) 


Up to eight (8) functions may be com- 
puted and plotted at one time. To plot 
more than eight functions, the fragment 
field file must be repositioned, and 
ancther set of plots requested. 


2) 


* ; 
Functions 08 and 09 must both be specified if! either is, and must 
be the last function codes specified. 


Note that a function code of 3j or 6j 
causes two plots to be generated, and 
that specification of functions 08 
and 09 counts for 3 functions, even 
though only two plots are generated. 


PLOT cause function values to be written to 
the plot file (file 2) for later 
processing by a contour plotting program. 


KWIKPLOT generate contour plots of functions 
values on the printer. 
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RANGE read plot range parameters from data 
card following this control card. The 
format of this card is given below. 


SOJAC read characters io be used for printer 
plotter (see KWIKPLOT) from data card 
following this control card. The format 
of this card is giver below. 


LIST list data read with this control card. 


ee 
a 


Data cards read after a $OUTPUT,TABLES statement. Between two and ss 
five cards are read for each table or formula coefficient group. 


COLS FORMAT NAME DESCRIPTION 


Heading Card a 


1 Tl IT type code. If IT=1,2 or 3, table IS is si 
to be read. If IT=4,5 or 6, coefficient 
group IS is to be read. 


IF IT=0 and 1S=0, the end of this group ag 
of data cards has been reached. 


2 Tl Is table or coefficient group number. a 
O<ISS9. 

3-5 (ignored) : 

6-10 15 NTB (used only for tables) number of M-V | 


pairs to be read. NTB<14. 


11-15 =o LL JDUM (used only for tables) interpolaticn 
mode. All table values are read in : 
true form, but linear interpolation of : 
the table during processing may require 
either a log or linear scale, according 
co JDUM: 
=00 M linear, Vv linear | 
=01 M lirear, V log | 
=10 M log, V linear 
=]11 M Jog, V log 


Coefficients (formula coefficients only) 


1-70 7E10.0 (COEF(1) ,I=1,7) - Thor formula coefficients, 
group 1s, 
Fragment Mass Values (tables only) 
NTE masses are read, 7 per card, 
1-70 7E10.0 (XTB(I) ,1=1,7) - Fragment mass values, (grains) 
for table IS. 
Damage Threshold Velocity Values (tables only) 
NTB velocities are read, 7 per card. 


1-70 7E10.0 (YTB(I) ,I=1,7) - Damage threshold velocities 
ft/sec) for table Is, 
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Data cards read after a S$OUTPUT,PARAMS statement. Between two and 
four cards are read fr each fragment field function. 


cOLs. 


FORMAT 


NAME 


DESCRIPTION 


Heading Card 


1 


Il 


Il 


Al 


IT 


Is 


IVP 


NF 


Function type. If IT=0 and 1S=0, the 
end of this group of data cards has been 
reached. If IT=0 and IS>0, then the 
function requires neither a table nor a 
Thor formula, and the function type is 
given by IS. IF IT20, then IT deter- 
mines the function type. 


=] number of damaging fragments per 
sq ft, using table IS. 

=2 probability of damage, using table IS. 

=4 number of damaging fragments per sq ft 
using coefficient set IS in the Thor 
formula. 

=5 probability of damage using coeffi- 
cient set IS in the Thor formula. 


table number er: IT=l or 2), coefficient 
group number (if IT=4 or 5), or function 
type (if IT#=0): 
=] totai number of fragments per sq ft 
=2 total fragment mass (1b) per sq ft 
=3 total fragment momentum (1lb-sec) 
per sq ft 
=4 total fragment kinetic energy (ft-Lb) 
per sq ft 


=8 average fragment impact velocity (ft/sec) 


=9 average fragment impact angle (radians) 
(ignored) 

Function variant code. ‘A blank is treated 
as a ‘zero’ character. IVP is used to dis- 
pace plots of the same function, but 
with different parameters, 

(ignored) 


Number of parameters (O<NF<14) If NFC14, 
parameters not read are set to zero, 
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COLS. FORMAT NAME DESCRIPTION 
Title Card 


1-72 12A6 (ITTLF (I) ,1=1,12) - function title, which will 
appear on all plots of this function, 


Parameter Card(s) 


NF parameters are read, 7 per card. If NF=0, no cards are read. 


1-70 7E10.0 (F(1) ,I=1,NF) - parameters affecting function 
performance. 
Parameters 1 through 7 have the same meaning for all functions: 


F(1) miminum fragment weight (grams). All lighter frag- 
ments are ignored, 


F(2) minimum fragment velocity (ft/sec). All slower 
fragments are ignored. 


F(3) minimum fragment momentum (lb-sec). All fragments 
whose momentum is too smail are ignored. 


F(4) minimum fragment kinetic energy (ft-lb). All frag- 
ments whose kinetic energy is too small are ignored. 


F(5) - F(7) (reserved) 


Parameters 8 through 14 are used only if IT»0. 
F(8) target plan area (£t*) 
F(9) average target elevation area (£t2) 
F(10) - F(11) are used only if IT=4 or 5. 
F(10) target thickness (in.) 
F(11) target hardness, where applicable. 
F(12) - F(14) (reserved) 
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coLs. 
1 


15 


16-20 
21-30 
31-40 


FORMAT 
Il 


12 


14 


Tl 


E10.0 
E10.0 


NAME 
IT 


Is 


ISM 


ISc 


OZMIN 
OZMAX 


Data cards read after a SOUTPUT,LIMITS statement. One card is read 
for each function. 


DESCRIPTION 


Function code (see PARAMS discussion 
above) 


subcode or table code (see PARAMS discus- 
sion above) 

Note that if IT=0 and IS=0 the end of 
this group of data cards has been reached. 
(unused) 


Plot smoothing code. IF ISM=0, do no 
smoothing. IF ISM=1 to 5, smooth by re- 
placing each value by the average of the 
values in a 2*ISM+l raster square cen- 
tered on that value, 


Plot scale code. ISC=0 to request a 
linear plot. ISC=1 to request a plot of 
the log of the function. 

(unused) 


minimum value to be plotted. 


maximum value to be plotted. 


If OZMIN>OZMAX, the plot range is set 
interna lly. 


coLs. FORMAT. NAME DESCRIPTION 


COLS. FORMAT. CNAME DRO CR ES 


Data card read after SOUTPUT,EXEC statement. See discussion of 
the EXEC parameter for the list of valid function codes. 


1-2 I2 first function code 
3-5 (ignored) 


Subsequent function codes go in columns 6-7, 11-12, etc. 


Data card read after SOUTPUT,SOJAC statement. 


1-5 (ignored) 
6-10 15 NZ number of function value ranges to be 
plotted on printer. NZ<50 
11-19 (ignored, leave blank) 
20-69 50Al (CH(1) ,1=1,NZ) - plot characters for each func- 


tion value range. 


Data card read after SOUTPUT,RANGE statement. 
1-5 (ignored) 


6-10 I5 NX number of resolution increments in X 
direction (tail to nose). NX<8l. 


11-15 15 NY number of resolution increments in Y 
direction (laterally from waist of muni- 
tion). Normally NY=(NX-1)/2+1. Nx<4l. 


16-20 (unused) 

21-30 £10.0 XMIN minimum X value (ft) 

31-40 £10.0 YMIN minimum Y value (ft). Normaliy YMIN=0. 
41-50 £10.0 DX X increment (ft). Normally NX,XMIN and ‘DX 


are set so that XMAX=-XMIN where XMAX 
is XMIN+DX*(NX-1). 


51-60 E10.0 DY Y increment (ft). Normally NY and DY 


are set so that YMAX=XMAX where YMAX 
is DY*(NY-1) if YMIN=0. 
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Sample Data Deck 
the deck listed on the next page is set up to do the foilowing 


job: 


1) Read munition data for munition No, 4 (mod 1). 
7 Gaps in the data are to be filled in, and the 
data is to be listed neatly in tabular form. 


2) A fragment field is to be computed, on a polar 
grid consisting of 8 ranges spaced 250 feet 
apart, and 18 azimuth rays, spread evenly over 
a half-circle (at azimuths 5, 15, ..-, 175 deg). 
After the fragment field has been written onto 
tape 4, an end of file is written and the tape 
is rewound. 


Baler er than RE Tie a a hee 


3) Output specifications are read. These include 
a. the printer-plotter character set, 
b. the size and resolution of the plot, 


ce. the trivial (O-velocity) fragment thresh- 
old velocity table, 


d. the function limits and scale (all are to 
be plotted in log units), and 


e. the function titles and parameters. 
4) Plots are generated. The function specifica- 


tions 01 and 30 cause functions 01, 10 and 20 
to be generated and plotted. 


j A-19 


SPRINT CONTROL+DATA 
SORNNANCE FORMATS2,NORM, SMOOTH, LIST 
MUNITION 94 MOD 4 

16 36 4-1 660, 


(The munition data goes here.) 


PROM NA LST ORNGEZOO. ONRNG SG $ 2000 FOOT RANIUS, 10 DEG aZ SECTORS 


GROUND LEVEL FRAGMENT &TELD 
TIVE 
STAPE 4,wE OF 
STAPE 4,26 NTN 
SOUTPUT SQUAG LIST 
9 


7 ~123495607094 
SOUTRPET RALGE LIST 
B10 44 -2u0c. he 50, bn, 
OUTPUT TABLES LIST 
10 2 g 
Qs inno, 
Ne Qa. 


EMD 
HOUTPUT KLImtts List Thre 
a1 4 100M16667 116667 [ReprodUCriie cop 
19 1 00016667 .16667 best ® 

20 1 NOL ft 


ae) 
BOUTPUT PARAMETERS LISI 


O4 4 
AREAL DENSITY OF ALL FIKAGMENTS (TARGET 1S SPHERE OF UNIT CROSS*SECTION) 


QO, 9. CG. 6 
40 9 
AREAL SENSITY OF NA“AGING FRAGMENTS (K.F, GT 58 FT=-LB) 
de oe oe 54, 
1.33 9. 
20 9 
PROPABILITY GF HIT BY DAMAGING FRAGMENT (K.E, GY 58 FT-LAD 
Qa, 0. 0. 58, 
14.33 9.0 

£ND 
SQUTPUT EXEC LIST KWIKPLOT 
GA 30 
TIME 
$SToOp 
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Sentinel Records 


The munition data file (if one is used), the fragment field file, 
and the output function file (if one is used) each consist of 
several groups of data. Each group is preceded by a sentinel 
record which identifies the group, as described in the following 
section, In addition a special sentinel record is written after 


the last data group. 


A sentinel record exists in two forms, one for formatted files 
(the BCD munition data file), and one for unformatted files (the 
binary fragment field file andthe binary output function file). 


Formatted sentinel record: 


char. 1-6 '“HEAD*' 


char. 7-14 an 8-character search key. (If fewer 
than 8 characters, blanks are added.) 


Search for a formatted data group using the statement; 
STAPE n,BCD ,SRCH=key , BKSP 
Unformatted sentinel record: 


word 1” 'eHEAD* ' 


word oie an 8-character search key. (If fewer 
than 8 characters, blanks are added.) 


Search for an unformatted data group using the statement: 


STAPE n,SRCH=key , BKSP 


Two words are needed when using a computer with four characters 


per word. 


we 
Two words are needed when using a computer with four or six 


characters per word. 
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Sentinel Record Keys 
1. A 3-character key of the form 
XXI 
is used for munition data. (BCD tape) 


XX = munition ID number. O0¢XX<99 
I =a modification code. Normally [=0. 


2. A 5-character key of the form 
XX LIK 
is used for fragment field tables. (Binary tape 4) 


XX,I as above 
J,K = fragment fie’d modification code to reflect 
changes in resolution, target plane altitude, 
barrier placement, etc.). Normally J=K=0. 
3. An 8-character key of the form 
XXIJKYYL 
is used for output function tables. (Binary tape 2) 


XX,I,J,K as above 
YY = output function code. (See $OUTPUT discussion.) 


L = output function modification code(to reflect 
changes in threshold or target parameters, etc.). 
Normally L=0, 


4. The sentinel record which should end a file has an 8-character 
key which is 
99999999 
A sentinel record may be written with the statement 


STAPE n{ ,BCD] ,WEOF 


A-22 


sesame meme ace See ca eb Ge BES ge OEE ss 


d 


7 
} 
Cad 


ERC Aina ore eee terremnnte ere wa  n  e  e 


Timing ; 
The generation of the fragment field is the most time consuming 
part of the program. Several parameters affect timing in this 
phase. 

NM number of mass categories into which munition 

data is divided 

NAZ number of azimuth angles 

NEL number of elevation angles 

DRNG range increment (ft) 


All operations are performed for each mass category at each azi- 
muth. Low-register fragment trajectories are computed for ranges 
R;=DRNG*i, out to R:=Riy> which corresponds to an initial eleva- 
tion angle of Onax' High-register fragment trajectories are com- 
puted for elevations a,=(i-1/2) 4a, where Aa=90°/NEL, and where 


A nays 44690 holds. 
Thus the computation time for low register fragment trajectories 
increases as DRNG decreases, and the computation time for high 


register fragment trajectories increases as NEL increases. So 
the expected computation time is given by 


T=A*NM*NAZ* (1+B/DRNG+C*NEL) 
Because all runs done for this project used the same values of NAZ, 


DRNG and NEL, the constants A, B and C cannot be determined at this 
time. However, the approximate relation (in sec) 


T=7 x NM 


and has been observed on the Univac 1108, with NAX=18, DRNG=250 
and NEL=18. 
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